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ABSTRACT 

Using the largest homogeneous quasar sample with high-quality optical spectra and robust radio 
morphology classifications assembled to date, we investigate relationships between radio and optical 
properties with unprecedented statistical power. The sample consists of 4714 radio quasars from 
the Faint Images of the Radio Sky at Twenty cm (FIRST) with ^20 > 2 mJy and with spectra from 
the Sloan Digital Sky Survey (SDSS). Radio morphology classes include core-only (core), core-lobe 
(lobe), core-jet (jet), lobe-core-lobe (triple), and double-lobe. Electronic tables of the quasar samples, 
along with spectral composites for individual morphology classes, are made available. We examine 
the optical colors of these subsamples and find that radio quasars with core emission unresolved by 
FIRST (on ~ 5" scale) have a redder color distribution than radio-quiet quasars (5'2o < 1 mJy); 
other classes of radio quasars have optical color distributions similar to the radio-quiet quasars. 
This analysis also suggests that optical colors of z < 2.7 SDSS quasars are not strongly (< 0.1 mag) 
biased blue. 

We show that the radio core-to-lobe flux density ratio (i?) and the radio-to-optical (?-band) ratio 
of the quasar core (i?i) are correlated, which supports the hypothesis that both parameters are 
indicative of line-of-sight orientation. We investigate spectral line equivalent widths as a function of 
R and including the [OIII] narrow line doublet and the CIV A1549 and Mgll A2799 broad lines. 
We find that the rest equivalent widths (EWs) of the broad lines correlate positively with i?i at the 
4(T-8(T level. However, we find no strong dependence of EW on R, in contrast to previously published 
results. A possible interpretation of these results is that EWs of quasar emission lines increase as the 
line-of-sight angle to the radio-jet axis decreases. These results are in stark contrast to commonly 
accepted orientation-based theories, which suggest that continuum emission should increase as the 
angle to the radio-jet axis decreases, resulting in smaller EW of emission lines (assumed isotropic). 
Finally, we observe the Baldwin effect in our sample, and find that it does not depend strongly on 
quasar radio morphology. 

Subject headings: quasars: general — quasars: emission lines 



1. INTRODUCTION 

The unification paradigm for radio-loud active galac- 
tic nuclei (AGN) suggests that many classes of 
observationally-distinct sources appear dissimilar be- 
• cause of orientation e ffects (e.g. Urry & Padovani 1995; 
' iJackson fc Walllfl999n . Because we cannot (yet) visit an 
I AGN or observe it from different angles, we often take 
a statistical approach to this problem, piecing together 
observations of many different sources to investigate this 
theory. Intrinsic source properties, such as age, size, and 
luminosity, must also play a role in the standard AGN 
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picture, complicating investigations of orientation mea- 
sures. The "holy grail" of evidence for the unification 
paradigm would be an observational parameter that is 
understood to correlate directly with orientation angle. 
Owing to its variety and complexity, quasar radio mor- 
phology is often considered to hold the key to unlocking 
the orientation mystery. 

Statistical studies of radio emission from extragalac- 
tic sources have recently entered a new era, thanks to 
the availability of large sky-area high-resolution radio 
continuum surveys that are sensitive to mJy flux den- 
sity levels, such as the Faint Images of the Radio Sky 
at Twenty cm fFIRST: iBecker et al.l [19951 ) survey. In 
this paper, we present a sample of 4714 radio quasars, 
spectroscopical ly-confirmed by the Sloan Digital Sky 
Survey fSDSS: iYork et al.ll2000l ) and with robust visual 
classifications of radio morphology from FIRST images. 
This is the largest sample of visually morphologically- 
classified radio quasars to date. We use the radio quasar 
sample to investigate the relationship between radio 
morphology and optical spectral parameters in the con- 
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text of orientation theories. 

In standard unification theory, the basic picture of a 
radio-loud quasar is that of a galaxy with an accret- 
ing central black hole (an AGN) surrounded by a dust 
torus in the plane of the accretion flow, and with power- 
ful relativistic jets directed outward along the rotation 
axis of the black hole. The anisotropic emission leading 
to orientation effects is thought to have two principle 
causes: obscuration of the central object and immedi- 
ate environs when viewed through the dust torus, and 
relativistic boosting of the jet emission when the line of 
sight is close to the jet axis. These effects have visible 
consequences in both the radio and the optical regimes. 
Simple models of AGN core emission motivate descrip- 
tions of how radio properties, such as morphology and 
spectral slope, depend on orientation. Theories of other 
orientation indicators, such as spectral lines, have been 
empirically motivated by observed correlations with the 
aforementioned radio parameters. 

Determining the orientation angle to any individ- 
ual AGN is extrem ely difficult, if not impossible. 
IGhisellini et al.l (|1993[ ) discuss a possible method, which 
requires high-resolution radio observations of the rela- 
tivistic jets. The viewing angle is related to the appar- 
ent jet velocity and the relativistic Doppler factor, S. An 
estimate of 6 is based on the "classical" condition that 
the synchrotron self-Compton flus0 from the jet can- 
not exceed the observed flux at X-ray frequencies. The 
self-Compton flux depends on S, on the angular size of 
the boosted core component, and on the Compton self- 
absorption frequency. The latter two properties can be 
determined from multi-frequency observations with high 
spatial resolution using Very Long Baseline Interferom- 
etry (VLBI); the necessary observations, however, are 
prohibitively expensive in terms of time and resources. 
Furthermore, Doppler factors (and therefore derived in- 
clination angles) measured on parsec scales can differ 
dramatically from those derived on kiloparsec scales ow- 
ing to even modest changes in jet acceleration or direc- 
tion (e.g., Harris & Krawczynski 2006). 

A common approach to orientation studies is to use 
large data samples and statistical indicators of ori- 
entation. One such parameter that has been em- 
ployed to great success is the radio core-to-lobe ra- 
tio, R. The extended emission from AGN radio 
lobes is thought to be isotropic, while emission from 
the relativistic jets is Doppler boosted. The param- 
eter R should therefore depend strongly on orienta- 
tion, wit h high R indicating a small angle to the 
iet axis (lOrr fc Browni 119821: iKapahi fc Saikial 119821: 
iMorisawa &: Takaharal 119871 ). Many other studies, 
considered both individually and collectively, support 
the use of i? as a sta t istica l indicator of orientation. 
lAntonucci &: UlvestadI ()1985f ) found that R correlates 
with optical polarization, vari ability of the co r e, and 
one-sided radio morphology. iZirbel fc BaumI (|1995f ) 
found that R correlates with total radio power in 
double-lobed radio sources, but with significant scat- 

^ "Synchrotron self-Compton" refers to inverse Compton scat- 
tering of synchrotron photons off of the electron population that 
is producing the synchrotron radiation. 



ter (2 orders of magnitude). iMorganti et al.l (|1997[ ) 
found higher values of R in sources thought to be 
close to face-on as deter mine d by radio morphol- 
ogv. iphiaberge et al.l (|1999f ) and iHardcastle fc Worralll 
(|200(F found that R correlates with core radio luminos- 
ity, while iKha rb fc Shastri (2004) found a correlation 
with core optical luminosity. 

The scatter in these correlations with R indicate that 
it is not a direct measure of orientation for any individ- 
ual source, and that some other factor or factors must 
also contribute. Age is a likely factor, as quasar radio 
lobes are expected to expand and d i m ov er time (e.g. 
iBlundeU et all fl99§. iBarai fc WiiTal [20071 ). Environ- 
ment may also affect extended radio emission, as demon- 
strated by the existenc e of HYbrid MOr phology Ra- 
dio Sources (HYMORS) (i Gopal-KrishmTfc Wiita 20Q3). 
HY MORS are radio-loud quasars having a different 
iFanaroff fc Rile y (1974) (FR) morphology on either side 
of the core. FR I sources have two visible jets that have 
lost their b ulk re lativistic velocities on large (> 10 kpc) 
scales (e.g.. iBick ncll 1995?) and that often terminate in 
diffuse lobes; by definition, the majority of their emis- 
sion arises from the inner half of their total extent. FR II 
sources are edge-brightened as their emission predomi- 
nantly comes from their lobes; they are thought to have 
jets that remain relativistic out t o very large distances 
and terminate in hot-spots (e.g., iGopal-Krishna et al.l 
12006( 1. FR I morphology is typical in sources with radio 
power P < lO^''^^^ WHz^^ while FR II sources typi- 
cally have P > 10^'*"^^ W Hz"^; the actual transition lu- 
minosity rises with optical luminosity of the host galax y 
(ILedlow fc Owed [1991 TGopal-Krishna fc Wiital [2001 . 
Environmental asymmetries could lead to different jet 
interactions on either side of the central source, even 
if the op positely-directed jets are similar in power and 
velocity ([Gopal-Krishna fc Wiital I2000D . The majority 
of quasars in the sample we present are, unfortunately, 
too distant for accurate determination of FR class from 
their radio images. 

It has been suggested that the radio-to-optical ratio 
of the quasar core, Ry, is a better statistical measure 
of core-boosting, an d theref ore of orientatio n, than R 
([Wills fc Brotherton 1995; K harb et al.|[20Tol) . i?v mea- 
sures core-boosting by using the core's optical flux to 
normalize the observed (boosted) core emission. In con- 
trast, R normalizes core-boosting using the extended 
radio emission; R is therefore sensitive to the jet's envi- 
ronmental interactions far from the nucleus, while Ry is 
not. It is important to note that using Ry as a statistical 
parameter assumes that the optical core emission origi- 
nates in the accretion disk. For face-on sources such as 
blazars, viz., optically-violent variables and BL Lac ob- 
jects, the beamed synchrotron emission dominates the 
optical spectrum. Such sources could therefore contam- 
inate our analysis. However, at such extreme angles, 
the beamed synchrotron should overpower the line emis- 
sion entirely; the fact that SDSS quasars are identified 
as having one or more very broad lines suggests that 
the optical emission is dominated by the accretion disk 
rather than a strongly beamed continuum. For this pa- 
per, we measure the radio-to-optical ratio of the quasar 
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core using the i band {'^ 7481 A). We will refer to this 
ratio as Ri rather than Ry. 

Orientation is thought to affect the optical spec- 
trtun throu gh orientation-dependent obscuration (e.g., 
INenkova et al.i |2008I . and references therein) and/or 
inclination-dependent emission from an accretion disk 
surrounding the central black hole of an AGN (e.g., 
lBakeii[l997l ). Evidence of these effects is typically found 
via optical parameters, such as emission line profiles or 
spectral index, that correlate with the known radio sig- 
nat ures of orientation: R and spectral slope. For exam- 
ple, iBaker fc HunsteadI (j 19951 ) made spectral compos- 
ites for sets of quasars grouped by R and found that as 
the supposed viewing angle to the jet axis increases, the 
optical continuum steepens, the Balmer decrement in- 
creases, and line widths increase. They concluded that 
reddening is strongest in quasars that are more lobe- 
dominated, in accordance with the basic unification the- 
ory 

In this paper, we revisit spectral line correlations with 
these two proposed orientation parameters, R and 
using a large quasar sample with well-determined core 
and lobe flux densities. The remainder of the paper is 
laid out as follows. In Section [21 we discuss the data 
selection and compilation of the radio quasars sample. 
In Section [31 we discuss how automated estimates of 
quasar morphology motivate the more detailed classi- 
fication presented in this paper, and describe the radio 
morphology visual classification method; we also investi- 
gate the intrinsic color distribution of quasars as a func- 
tion of radio morphology. Section [H investigates how R 
and Ri are correlated. Section[5lpresents a set of quasar 
spectral composites based on morphology, on R, and on 
Ri. In Section [6l spectral line properties in individ- 
ual spectra are investigated quantitatively as a function 
of the two proposed orientation parameters. We sum- 
marize our results in Section [71 Throughout, we use a 



the sky at 20 cm (1.4 GHz) with a beam size of 5'.'4 and 
an rms sensitivity of about 0.15 mJy beam^^. FIRST 
observed 9,000deg^ at the north Galactic cap and a 
smaller ~ 2.5° wide strip along the Celestial Equator 
in the South Galactic cap, overlapping the SDSS foot- 
print. FIRST is 95% complete to 2 mJy and 80% coin - 
plete to the survey limit of ImJy ([Becker et al.lll995l ). 
The survey contains over 800,000 unique sources, with 
astrometric uncertainty of < 1". From fitting a two- 
dimensional Gaussian to each co-added image, FIRST 
measured the peak flux density at 20 cm (S'pcak) and the 
integrated flux density at 20 cm (6*20) for each source. 

2.2. The Sloan Digital Sky Survey quasar catalog 

We draw our quas ar sample from the quasar 
catalog presented by [Schneider et al.l ([2007[ . here- 
after S07), which is based on the fifth data re- 
lea^ (DR5) of the SDS^ (lAdelman-McCarthv et all 
[2007f) . Tec hnical det ai ls ab o ut the surv e y ca n 
be found in lYork et all ([200l: IGunn et al.l (l2006l). 



Fukugi ta etal.1 ([19961) : [ Gunn et al.) ([1998D: [Hogg et all 



20011: [Lupton et al.[ ([2002D: [Scranton et al.l ([200li 



iSmith et al. (20021: Ivczic et al.' j2004| ): [Tucker et all 

m 



cosmology with Hq 
Qa = 0.7. 



2. DATA 



([2006 '): Padmanabhan et al. ( 2 00a) ■ The astro metric 
accuracy of the SDSS is < 0.1" ([Pier et al.ll200l . 

A small fraction of SDSS photometric sources are 
selected for spectroscopic observation according to 
several targeting algorithms. The quasar algorithm 
(jlichards et al. 2002'!) targets afl 15 < i < 19.1 
(for «-band apparent magnitude i) unresolved sources 
within 2" of a FIRST radio detection. Some quasars 
are targeted based on their broadband SDSS colors; 
others arc targete d fortuitously via the algorithms 
for gala xies ( Eisens tein et al.l [20011 : iStrauss et al.l 120021 : 
iBlanton et al.ll2003l) . 

The S07 sample contains 77,429 sources from w 
70kms"^ Mpc"\ f^A/ = 0.3, and 5740deg^ of sky with i-band absolute magnitude Mi < 

—22 and at least one emission line having FWHM larger 
than 1000 km s^^, or which contain complex, broad 
absorption features. Rest-frame i-band absolute mag- 
nitudes were calculated assuming a power-law spectral 
energy distribution with spectral index ctopt = —0.5. 
An updated quasar cat alog has recently been published 
([Schneider et al.l I2010L here after SIO) from the S DSS 
seventh data release (DR7; lAbazaiian eFaDHQOi). In 
the updated DR7 catalog, 171 objects from DR5 were 
dropped, typically because of small changes in photom- 
etry or because their spectra do not satisfy the mini- 
mum line wid th as measured using a principal compo- 
nent analysis ([Schneider et al.[[2010l ). We exclude these 
171 objects from our analysis. We use the redshifts, 
astrometry, and photometric measurements from SIO. 

2.2.1. Spectral line measurements 

To obtain spectral line measurements, we use the 
catalog of SDSS DR7 quasar properties presented by 
iShen et all ([20101 ). The catalog provides emission line 
measurements for Ha, H/3, CIV A1549, Mgll A2799, 



We describe below the radio and optical sky surveys 
from which we draw our quasar sample. The radio sur- 
vey, FIRST, is the deepest large-area radio sky sur- 
vey undertaken to date. Performed at 20 cm, it is 
sensitive to the high-frequency quasar core emission; 
the survey contains high-resolution radio images. The 
SDSS provides the largest-ever homogenous sample of 
spectroscopically-confirmed quasars. FIRST was de- 
signed to cover the same region of sky as the SDSS, 
making this combination of sky surveys a powerful tool 
for statistical studies of radio quasars. 

2.1. FIRST 

The FIRSTCI survey ([Becker et all [Ml) used the 
NRAC|3 Very Large Array (VLA) telescope to observe 



^ Ihttp: //sundog. stsci . edu| 

^ The National Radio Astronomy Observatory is a facility of the 
National Science Foundation operated under cooperative agree- 
ment by Associated Universities, Inc. 



' |http: //www. sdss ■ org| 
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and the [O III] doublet at A4959 and A5007 (among oth- 
ers). Each line was modeled using one or more Gaussian 
profiles, as appropriate; broad and narrow components 
of the emission lines were fit separately. This method 
is more robust than that used by the SDSS pipeline, 
which fits a single Gaussian to each line profile. The 
Shen et al. fits deteriorate for low-quality spectra; for 
the present paper, we limit our emission line analysis to 
spectra with median signal-to-noise ratio per pixel (69 
km s-i SDSS pixels) > 5. 

The presence of broad absorption lines (BALs), 
which are evidence of quasar o utflows and are at leas t 
2000 km s-i wide by definition (iWevmami et al.lll99lD. 
can s kew spectral line measurements. I Gibson et al.l 
(|2009( ) visually examined each DR5 SDSS quasar spec- 
trum and identified those with BALs near the Si IV, 
CIV, Al III, and Mgll lines. For analysis of individual 
lines, we exclude quasars with a BAL associated with 
that wavelength. Because the presence of a BAL can 
also affect the spectral shape in the ultraviolet regime, 
we exclude every BAL quasar from our analysis of spec- 
tral slopes in Section 13.41 we also exclude BALs from 
our composite spectra in Section [5. II 

3. RADIO MORPHOLOGY OF QUASARS 

Radio morphology is crucial for categorizing quasars 
and understanding their underlying physical properties. 
The apparent morphology of a radio source is deter- 
mined by intrinsic parameters such as age, size, and jet 
power, for example, but also by orientation. In this 
section, we present the classification of radio quasars 
by radio morphology and define the morphology classes 
used throughout the rest of this paper. We then discuss 
the fraction of highly-reddened quasars as a function of 
morphology class. 

3.1. Lessons from automated morphology classification 

We begin by giving a brief overview of work by 
iKimball fc Ivezid (|2008i . hereafter KI08) to classify ra- 
dio sources by morphology and to describe the optical 
and radio properties of the different morphology classes. 
Their work motivates the more detailed morphology in- 
vestigations we now pursue. 

KI08 compiled a large, homogeneous catalog of ra- 
dio and optical sources for the purpose of exploring 
the multi-wavelength properties of AGNs. A simple 
method was used for automatic morphology classifica- 
tion based on flux densities obtained at varying spatial 
resolution from different radio sky surveys. They used 
this method to deflne three morphology classes: "com- 
pact" (unresolved at 5"), "resolved" (resolved at 5" but 
fairly concentrated), and "complex" (having signiflcant 
extended emission beyond 5"). KI08 found that these 
three classes, defined entirely by their radio emission, 
correspond to different types of optical sources (with 
significant overlap among the classes). For example, 
the complex and resolved classes consist of optically- 
identified galaxies with steep radio spectra, while flat- 
spectrum quasars typically fall into the compact class. 
These results are consistent with the basic unification 
theory of radio-loud objects, which suggests that radio 
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Fig. 1. — Examples of the different morphology classes as labeled 
(see text for definitions) in FIRST 2' X 2' stamps normalized to the 
brightest pixel and scaled linearly. FIRST instrumental resolution 
is ~ 5". The fourth and fifth rows show sources with increasing 
core-to-lobe flux density ratio from left to right. 

AGNs viewed along the jet axis appear as compact and 
flat-spectrum in the radio (because core-boosting dom- 
inates the radio emission) and as quasars in the optical 
(because the central object is unobscured). Conversely, 
AGNs viewed perpendicular to the jet axis tend to be 
dominated by steep-spectrum extended radio emission, 
and appear as galaxies in the optical regime. 

With the new quasar classifications presented here, 
we extend the analysis of KI08 by examining the depen- 
dence of quasar spectral properties on radio morphology. 
Additionally, our classifications separate the "complex" 
class of KI08 into finer morphology categories. 

3.2. Radio morphology classification of quasars 

3.2.1. Sample Definitton 

We defined our initial quasar sample using 77,258 
quasars from S07 (the full catalog of 77,429 sources 
from the SDSS DR5 catalog, less the 171 sources tha t 
were dropped from the updated DR7 catalog; see ^2.2\i . 
We divided the S07 sample into a main sample of "ra- 
dio quasars" and a secondary sample of "radio-quiet 
quasars". The radio sample consists of 4714 quasars 
with a nearby FIRST source whose integrated flux 
density is 5*20 > 2 mJy, in order to ensure a suffi- 
ciently bright sample for reliable morphology classifica- 
tion. The radio-quiet sample comprises 65,253 quasars 
that lie in the area of FIRST coverage but do not have 
a FIRST counterpart H We exclude from our analysis 

The term "radio-quiet" has more than one definition in the 
scientific literature. It has been used in reference to a limiting 
radio luminosity (e.g., Dun lop et al. 1993), a s well as t o a limiting 
radio-to-optical flux density ratio (e.g., Peaco ck et al.| [T986). The 
definition employed here is different from both of these, as our 
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an additional 5525 quasars that lie outside the FIRST 
survey area. We also exclude 1449 "faint radio" quasars 
with a faint {S20 < 2 mJy) FIRST counterpart within 
2". 

The radio sample may contain a small number of 
false matches. The distribution of FIRST— SDSS po- 
sitional offsets is a Rayleigh distribution with a = 0.2 
and a small tail out to 2" ([Schneider et al.l [20070 . in- 
dicating that the number of spurious FIRST — SDSS 
matches within 2" is small. The FIRST survey den- 
sity is ~ 63 deg^^ for sources brighter than 2 mJy. It 
follows that the false matching rate in the S07 catalog is 
6.1 X 10~^ per source. We therefore expect no more than 
4-5 spurious radio sources in our radio quasar sample. 

3.2.2. Classification of radio images 

We classified the main sample by visual examination 
of 2' X 2' FIRST images centered on the quasar opti- 
cal position. The great majority of radio quasars are 
smaller than 2' in angular extent Me Vries et al.l [20061 
hereafter dV06). For images that showed emission > 1' 
from the optical position, we also reviewed a 4' x 4' field. 
We categorized the images into five morphology classes: 

1. core, a quasar with radio emission only at the op- 
tical position; 

2. jet, a source with radio emission from the quasar 
core and a radio jet; 

3. lobe, a source with radio emission from the quasar 
core and from a lobe; 

4. triple, a quasar with lobe-core-lobe morphology; 

5. knotty, a single source with many emission regions 
or a complex emission map. 

Example images of the first four of these classes are 
shown in Figure [Tj 

Two of us (AEK and PJW) examined the 4714 im- 
ages, and agreed on the morphology classification of 
4460 sources (95%). We assigned these sources a qual- 
ity flag value of 1. Author ZI examined the remaining 
254 sources and acted as a tie-breaker in 175 cases; we 
assigned these sources a quality flag value of 2. The 
remaining 84 images were sufficiently complicated that 
either all three examiners classified them differently, or 
agreed that a definite classification is not possible; these 
objects are listed as "unclassified" with a quality fiag 
value of 3. The number of objects in each morphology 
class is listed in Table [1] The main radio quasar sample, 
including morphology classifications, is given in Table [2| 

3.2.3. Redshift Selection Effects 

The radio quasars range from 0.78 to 5.4 in redshift. 
Redshift can affect the morphology classification be- 
cause faint radio components drop below the FIRST de- 
tection limit at large recessional distance. Furthermore, 

"radio-quiet" sample is defined by falling below the limiting radio 
flux density of FIRST: ~ 1 mJy. 



lAriljtlj i 




Quasar Classes 




Sample 


Sample size 


All quasars 


77,258 


Excluded: outside FIRST area 


5525 


Excluded: 1 mJy < 520 < 2 mJy 


1449 


Radio- quiet 


65,253 


Double-lobed w/o a radio core 


317 


Main radio sample 


4714 


Core 


3433 


Unresolved 


2802 


Resolved 


631 


Jet 


183 


Lobe 


387 


Core-dominated 


233 


Lobe-dominated 


150 


(undetermined) 


4 


Triple 


619 


Core-dominated 


164 


Lobe-dominated 


226 


Irregular 


129 


(undetermined) 


100 


Knotty 


8 


unclassified 


8 



extended source components may be unresolved at high 
redshift. Figure [2] shows the redshift distribution for 
each morphology class, and verifies that the fraction of 
sources with extended morphology (lobe, jet, or triple) 
is higher at low redshifts. 

Core sources can be divided into two categories, based 
on whether a source is resolved or unresolved in FIRST. 
From KI08, we define a dimensionless concentration pa- 
rameter as the ratio of integrated flux density (6*20) to 
peak flux density (Speak) at 20 cm, according to the 
formula 



"Resolved" core sources have log(0^) > 0.05 and "unre- 
solved" sources have log(0^) < 0.05. The resolved sam- 
ple likely consists of quasars with extended radio emis- 
sion that are still young, and therefore small, and/or at 
high-redshift such that the two radio lobes appear as 
one component. The unresolved sample should contain 
quasars with boosted radio cores, which will appear as 
point sources, and that do not have detectable extended 
emission. However, some of the unresolved sources, like 
the resolved cores, may instead be distant quasars whose 
extended emission is not resolved on the ^5" scale. 

3.2.4. Sub- classification of Triple and Lobe Sources 

To measure the core-to-lobe ratio R, we must deter- 
mine the flux densities of individual core and lobe com- 
ponents in the triple and lobe morphology classes. We 
identified these components directly from the FIRST 
images, where possible. Positive identifications were 
possible in 519 (383) out of 619 (387) sources with triple 
(lobe) morphology. In the remaining 100 (4) images, the 
field is too crowded to positively identify the lobe com- 
ponents associated with the quasar. For sources with 
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multiple emission points in a lobe (162 triples; 22 lobes) 
we use the sum of contributing components. 

We define "lobe-dominated" sources (lobes brighter 
than the core), "core-dominated" sources (lobes fainter 
than the core), and "irregular" sources (one lobe 
brighter than the core and the other fainter; this cat- 
egory only applies to triple-morphology sources). The 
bottom three rows of Figure [T] show examples of these 
three sub-classes among the triples. The number in each 
subclass is listed in Table [TJ The 20 cm flux densities of 
the cores and lobes are listed in Table [2] 

3.3. Douhle-lohed quasars without a radio- detected core 

Not all quasars with extended radio emission have an 
observed radio core. Unification theory suggests that 
such sources have jets aligned close to the plane of the 
sky. This population is excluded from our main ra- 
dio sample by the requirement of a FIRST counterpart 
within 2" of the quasar's op tical posit i on. T he fraction 
of omitted sources is small: iLu et al.l ()2007[ ) estimated 
that 8% of radio quasars are missed when using a 2" 
optical — radio matching radius. These so-called douhle- 
lohed quasars are useful to test for selection bias in the 
SDSS spectroscopic targeting algorithm fi i3.4p . 

To find a sample of double-lobed quasars, we started 
with our radio-quiet sample (quasars having no FIRST 
counterpart within 2"). Potential lobe components were 
FIRST sources with no optical counterpart within 3", 
and which were not associated with one of the quasars 
in our triple or lobe classes. We matched the radio-quiet 
quasars to these FIRST sources within 120"(dV06). For 
simplicity, we avoided crowded fields by limiting our fo- 
cus to quasars with no more than two potential FIRST 
lobes within 120". Among the remaining candidates, 
we required a lobe-corc-lobe opening angle of greater 
than 120° and a lobe- lobe flux density ratio between 0.1 
and 10. These choices of constraints were motivated by 
the selection algorithm for optically-faint triples (fJU de- 
scribed in detail in the Appendix). Our criteria result 
in a sample of 317 double-lobed quasars. 

A different, more extensive method to identify double- 
lobed quasars even in very crowded environments is de- 
scribed by dV06. For simplicity, we limit our selection 
to quasars in sparse environments, such that the lobes 
are easy to identify. (dV06 identified optical positions 
of their candidates, but not the individual lobe compo- 
nents.) We compare our sample to the dV06 double- 
lobed quasar catalog to estimate our completeness. The 
dV06 sa mple was drawn from the SDSS Third Data Re- 
lease (|Abazaiian et al.|[2005D and includes 143 quasars 
that are also in the S07 catalog. Our double-lobed sam- 
ple and the dV06 sample have 66 objects in common. 
The remaining 251 quasars in our sample are from later 
SDSS data releases than the dV06 sample. Of the 77 
sources which are in dV06 but not in our sample, 12 
were excluded because one of the potential FIRST lobes 
was associated with either an optical source or the ex- 
tended emission of a triple-morphology quasar. The re- 
maining 63 dV06 quasars were excluded from our sam- 
ple because they lie in a crowded field (more than 2 
FIRST matches within 120"). We successfully recov- 
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Fig. 2. — Redshift distributions for the radio quasars according 
to morphology. Left panel: Sources in each class as a function of 
redshift. Right panel: Ratio of resolved core, jet, lobe, and triple 
sources to unresolved core sources as a function of redshift. 

ered 66 dV06 quasars and miss no more than 63; the 
lower limit on the completeness of our sample is approxi- 
mately 50%. We note that we are biased against quasars 
in crowded radio sky regions, and those with lobes that 
were resolved into separ ate FIRST detec tions. 

We use the results of ILu et all ()2007D as a separate 
estimate of the completeness of the double-lobe sam- 
ple. Lu et al. found that 8% of radio quasars do not 
have a radio-core within 2" of the optical position. The 
size of our main sample is 4714 quasars, suggesting that 
there are ^ 410 corresponding FIRST — SDSS quasars 
that have only extended radio emission. We found 317 
double-lobed quasars, implying that our selection cri- 
teria have excluded a further ^ 93 sources. Based on 
these values, our double-lobed sample completeness ap- 
proaches 77%. These sources are listed in Table [H 

3.4. Fraction of highly-reddened quasars as a function 
of radio morphology 

Previous observations of quasar colors suggest that 
radio qu asars have redder S PSS colors than radio-quiet 
quasars (|Ivezic et al.l[2()()l . SDSS quasar targets are 
identified by selec ting outliers from the stellar locus in 
SDSS color space [Richards et"aLl (|2002[ ). Quasars with 
significant dust-reddening may lie on the stellar locus, 
and thus are not targeted by the color criterion. How- 
ever, every point source brighter than i = 19.1 is tar- 
geted if it has a FIRST counterpart within 2". Thus, it 
has been suspected that the difference in colors between 
radio and non-rad io quasars is merely a selection bias 
(llvezic et al.ll2Q0l . 

With our morphology classifications in hand, we 
can now examine these color differences in more de- 
tail, and quantify the optical color biases in the SDSS 
quasar sample. For example, systematic color differ- 
ences among the radio morphology classes would indi- 
cate a physical cause, because all were targeted with 
the same algorithm. However, if some radio quasar 
classes have the same color distribution as the radio- 
quiet sources, then it is unlikely that the colors of the 
radio-quiet quasars are strongly biased blue. 

We compare quasar ultraviolet spectral index (a^) 
distributions, as de termined with the method of 
[Richards et al.l ()2003[ ). who showed that the photomet- 
ric spectral index of quasars correlates strongly with the 
photometric 5 — i color. They defined the A ((7 — i) "color 
excess" as the difference between a quasar's g — i and 
the median quasar g — i at the same redshift. The color 
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TABLE 2 

Quasar, sample with radio morphology classifications 



SDSS name 


R.A.'' (J2000) Dec'' 


z 


class^ 


subclass*^ 


flag'^ 


'5'corc 


'^lobol'* 


•S'lobc2'^ 


000017.38—085123.7 


0.072423 


-8.856608 


1.2491 


J 




1 


9.77 






000028.82-102755.7 


0.120087 


-10.465497 


1.1518 


C 


U 


1 


2.23 






000050.60-102155.9 


0.210837 


-10.365531 


2.6404 


C 


U 


1 


20.39 






000051.56+001202.5 


0.214855 


0.200707 


3.9713 


C 


U 


1 


2.99 






000054.96+010143.4 


0.229009 


1.028724 


1.4617 


c 


u 


1 


3.74 






000111.19-002011.5 


0.296662 


-0.336539 


0.5179 


J 




1 


25.00 






000221.11+002149.3 


0.587988 


0.363706 


3.0699 


c 


u 


1 


13.88 






000222.47-000443.5 


0.593646 


-0.078752 


0.8106 


c 


u 


1 


3.89 






000442.18+000023.3 


1.175791 


0.006480 


1.0068 


L 


(J 


1 


3.71 


3.45 




000507.05-101008.7 


1.279415 


-10.169089 


1.2953 


C 


R 


1 


94.62 






000608.04-010700.7 


1.533519 


-1.116869 


0.9486 


T 


L 


1 


4.06 


51.14 


29.79 


000622.60-000424.4 


1.594198 


-0.073455 


1.0377 


C 


R 


1 


3879.24 







Note. — This subset of the table demonstrates format and content. The complete table is available in the electronic 
version of this paper. Flux densities are given in mjy units. 

SDSS position in units of decimal degrees. 
^ C: core; J: jet; K: knotty; L: lobe; T: triple; X: unclassified; — : no subclass. (Section 13.21 1 

U: unresolved core sources; R: resolved core sources; C: core-dominated triple or lobe sources; L: lobe-dominated 
triple or lobe sources; I: irregular triple sources; X: undetermined subclass for triples/lobes; — : source class has no 
subclass 

^ A flag of 1 indicates a source with recognizable morphology, such that the two first-round viewers initially agreed on 
the category. A flag of 2 indicates a more difficult case, where a third viewer's opinion was necessary as a tie-breaker. 
A flag of 3 indicates a source with undetermined morphology. 

° Total observed 20 cm flux density of the extended emission (lobes) of triple- and (o6e-morphology sources. Sources 
with undetermined lobe components (due to crowded FIRST images) are given a default value of -99. 



TABLE 3 

317 Candidate double-lobed quasars without a detected radio core 



SDSS position 



FIRST matches 



SDSS name 


R.A. ( 


;j2000) Dec 


z 


R.A. (J2000) Dec 


Score [mJy] 


000657.63-010358.8 


1.74016 


-1.06632 


1.4358 


1.73900 


-1.06724 


16.17 










1.74125 


-1.06559 


22.80 


001138.43-104458.2 


2.91015 


-10.74952 


1.2735 


2.91197 


-10.74503 


27.24 










2.90913 


-10.75229 


60.72 


002909.13-092110.1 


7.28805 


-9.35281 


1.5390 


7.26567 


-9.34251 


1.88 










7.30512 


-9.36770 


9.34 


002948.54+004447.5 


7.45228 


0.74655 


1.0078 


7.45183 


0.74489 


80.55 










7.45288 


0.75012 


36.91 


003236.76-001446.6 


8.15318 


-0.24629 


1.9073 


8.14862 


-0.23908 


8.15 










8.17664 


-0.26754 


1.17 



Note. — The full version of this table is available in the electronic version of this paper. Right 
ascension and declination are given in units of decimal degrees. 



TABLE 4 
Results of KS tests comparing 

aj/^ distributions of radio 
quasars and radio-quiet quasars 



Radio morphology 



P-value 



unresolved core 
resolved core 
jet 
lobe 
triple 
double-lobe 



3.8e-35 
0.36 
0.52 
0.54 
0.078 

0.0051 



Note. — P-values indicate the prob- 
ability that the radio quasar class has 
the same parent distribution as the 
radio-quiet class. 



excess correlates with ultraviolet spectral index, a]j 
as 



uv 



a„ 



= -0.6 - 2 X A{g - i) 



(2) 



(G. Richards, private communication); the correlation 
is not a function of redshift. The median scatter in a]l^ 
about the correlation is approximately 0.1. 

The top panel in Figure [3] compares the distribu- 
tions of radio-quiet quasars and radio quasa r s brig hter 
than i = 19.1, reproducing the llvezic et al.l (|2002[ ) re- 
sult with a much larger sample. The long tails to- 
ward negative values indicate sources with extremely 
red spectra. The radio quasars are shifted toward red- 
der colors: the radio-quiet sample has a median value of 
= —0.67 and the radio sample has a median value 
of = -0.75. 

The lower panel of Figure [3] shows cumulative dis- 
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Fig. 3. — Distributions of ultraviolet spectral index aj/^; more 
negative values indicate a redder continuum. Top; Normalized 
histograms of aj/^ for radio-quiet quasars (orange) and radio 
quasars (black) for sources with i < 19.1. Spectral index was 
calculated from g — i color-excess (see text). Dashed lines show 
the median values for each sample. The error bar in the lower 
left shows the median scatter of about the (5(g — i) versus 

oj/^ correlation. Bottom: Cumulative distributions for individual 
radio morphology classes and the radio-quiet sample. 

tributions of aj/^ for the individual radio morphology 
classes, including the double-lobed quasars. For sources 
that have core radio emission, vife apply a limit in radio 
flux density of 9. 1 mJy (14.0 mag, u sing the AB mag- 
nitude system of lOke fc GunnI [l983f) : this is the limit 
above which the resolved/unresolved separation for the 
core sources is most effective (K108). The main result of 
this analysis is that unresolved core quasars have a sig- 
nificantly higher fraction of extremely red sources than 
all other quasar classes. For example, unresolved cores 
make up 80% of the radio sample with aj/^ < —1.5, but 
only 45% of the radio sample with aj/^ > —0.5. 

To quantify the statistical significance of differences 
in Q^'^ distributions for radio quasars and radio-quiet 
quasars, we employ the Kolmogorov-Smirnov (KS) test. 
The P-value resulting from such a test represents the 
probability of obtaining the observed distributions by 
chance, if the two samples share the same parent dis- 
tribution. Table |4] shows the results of KS tests com- 
paring the radio-quiet quasar distribution to each radio 
class. The KS tests quantify what is shown qualita- 
tively in Figure [3l the only quasar class with an ex- 
tremely different [P ^ 0.01) color distribution from the 
radio-quiet quasars is the unresolved core class. This 
difference must be a physical effect because unresolved 
core sources were targeted with the same algorithm as 



resolved cores, triples, lobes, and jets. Because unre- 
solved core quasars dominate the radio sample, a com- 
parison of color distributions for the radio quasars and 
the radio-quiet qu asars shows the difference reported by 
llvezic et al.l (j2002[ ). Double-lobed quasars were targeted 
using the same criteria as radio-quiet quasars. Because 
those two classes have not dissimilar observed color dis- 
tributions, it is possible that they were drawn from sim- 
ilar parent color distributions. 

The simplest interpretation of our results is: (1) un- 
resolved core quasars are the only subsample of radio 
quasars that has an intrinsically redder color distribu- 
tion than other quasars, and (2) optical colors of SDSS 
quasars are not strongly (< 0.1 mag) biased blue. We 
note that these results are specific to quasars in the red- 
shift range z < 2.7; quasars in Figure[3]are roughly lim- 
ited to that range because of the redshift distribution of 
the S07 quasar catalog. 

It is not obvious why the unresolved cores have a 
much higher fraction of extremely red sources than the 
other samples. Figure [2] suggests that this result is not 
a redshift effect, because lobe and jet sources have simi- 
lar redshift distributions to the unresolved core sample, 
while triples and resolved cores have simi lar redshift dis- 
tribut ions to each other. As discussed in [Richards et al.l 
(|2003f) . dust-reddened spectra are not easily character- 
ized by a single power law. Our current analysis does 
not distinguish between spectra that are intrinsically red 
and those that are dust-reddened. We do not expect 
that the highly reddened sources in the unresolved core 
sample are caused by dust reddening, because these are 
sources that should have the smallest viewing angles and 
thus the least obscuration by the AGN torus. 

4. COMPARISON BETWEEN TWO CORE-BOOSTING 
PARAMETERS, R AND i?i 

Recent arguments suggest that R\ is a stronger sta- 
tistical measure of core-boosting, and therefore of orien- 
tation, than i? (see discussion in iJT]). i? normalizes the 
core boosting using the extended radio emission, which 
depends on age and environment in addition to orien- 
tation. These factors should have less influence on 
which uses the core's optical luminosity to normalize the 
core boosting. While a beamed synchrotron component 
(such as in blazars and BL Lac objects) would contami- 
nate the R\ estimates in face-on sources, the presence of 
broad lines in the spectra (required for inclusion in the 
SDSS quasar sample) suggests that such a component 
is not seen in our sources. In other words, it is unlikely 
that we have any directly face-on sources, as such ori- 
entations should result in an absence of any emission 
lines. 

We now compare values of R and R\ measured from 
our sample. If the two parameters were both perfect 
measures of orientation, then they would correlate per- 
fectly. However, factors that affect R and Ri separately, 
such as environment or age, can induce scatter in the 
correlation. 

4.1. Core-boosting parameters 

We define R as the core-to-lobe flux density ratio at 
20 cm, converting to the rest frame according to the 
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Fig. 4. — Correlation of R and Rj for 519 triple quasars with 
well-determined core and lobe flux densities. Panel A: Symbols are 
color-coded according to 20 cm core flux density. A magnitude of 
15 corresponds to 3.63 mjy; a magnitude of 10 corresponds to 363 
mjy. Panel B; R distributions of optically-bright and optically- 
faint triples, limited to sources with single-component FIRST de- 
tection of each lobe (see text). Solid line shows the observed-frame 
R distribution of points in panel A; dotted line shows observed- 
frame R for optically-faint triples. Panel C: Ri distribution of 
optically-bright triples (blue asterisks) and of the core and lobe 
sources (orange diamonds). 



K-correction formula 



R = R 



Sr, 



[20 cm] 



Si 



obc 



-(1 + ^) 



-aiobo) 



(3) 



where S'corc and S'lobo refer to observed 20 cm flux den- 
sities. We assume a core spectral index of Ofcorc = and 
a lobe spectral index of aiobo = —0.8 (KI08). 

The radio-to-optical ratio is typically referred to in the 
literature as Ry, the ratio of radio to visible light. For 
our purposes, we use the ratio of core 20 cm luminosity 
to the luminosity in the i band, which plays a large 
role in SDSS quasar selection. Formally, the «-band, 
with an effective wavelength of 7481 A, is slightly redder 
than what is typically considered the long-wavelength 
boundary of the optical regime at 7000A. We refer to 
this ratio as defined as 



log(i?i) = log 



L 



radii 



^optical ) 



M, 



radi( 



M,, 



-2.5 



(4) 



where M,; is the Galactic reddening corrected 
(jSchlegel et al.l 119981 ) and K-corrected i-band absolute 
magnitude from S07, who adopted a spectral index of 
ttopt = —0.5 to perform the K-correction. Mradio is 
the analogous "radio absolute magnitude" at 20 cm, ob- 
tained according to the formula 

Afradio = "^corc - k - D, (5) 

where D is t he distance modu lus and k represents the 
K-correction (|Hogg et al.ll2002[ ): 



k = -2.5 X (1 + acorc) X Iog(l -t- z). 



(6) 



The "radio apparent magnitude" Wcore is the observed 
20 cm flux density converted to a ppare nt magnitude on 
the AB^ system of lOke fc Grn^ (|1983D . 



m-corc = -2.5 log 

As before, we adopt acorc = 



3631 Jy 



(7) 



4.2. 



Evidence of a physical correlation between 
R and Ri 

Figure S] shows that R and i?i are strongly correlated. 
It also demonstrates that quasars with faint observed 
radio cores (indicated by symbol colors) have low core- 
boosting parameters, and vice-versa. The relationship 
between core brightness and core-boosting parameter is 
most likely a selection effect. To be classified as a triple, 
both of a quasar's lobes must be above the FIRST de- 
tection limit. Therefore, triples with high R must have 
cores that are much brighter than the FIRST limit. 

We now investigate whether the R — i?i correlation is 
a selection effect or an intrinsic property of quasars. 
For example, selection biases against high-i?i — low-i? 
sources (lower-right corner of panel A) or against low- 
Ri — high-i? sources (upper-left corner of panel A) would 
induce a spurious observed correlation. The evidence 
presented in panels B and C, discussed below, suggests 
that the correlation is a true physical effect. 

Panel B indicates that the dearth of sources in the 
lower-right corner of panel A is due to intrinsic quasar 
properties, for the following reasons. At a given ra- 
dio core fiux density, sources with high Ri are more 
likely to fall below the faint limit for SDSS spectroscopic 
targeting; we refer to these sources as "optically-faint 
triples", in contrast to "optically-bright triples" in our 
main sample. Optically-faint triples having lower R val- 
ues than optically-bright triples would be an indication 
that the empty lower-right corner in panel A is due to 
a selection bias. Panel B compares the R distributions 
of optically-bright and optically-faint triples. The cre- 
ation of the latter sample is described in detail in the 
Appendix. For simplicity of sample selection, we require 
the optically-faint triples to have each of their two lobes 
detected as a single component in FIRST. Therefore, 
in order to compare two equivalent data sets, panel B 
is limited to single-component lobes for the optically- 
bright triples as well. For the optically-faint triples, we 
report R values in the observed frame because we do 
not know their redshifts. Therefore, panel B also in- 
cludes the observed-frame R distribution for optically- 
bright triples. The distributions indicate that optically- 
faint triples have R values just as high (actually, even 
higher) than the optically-bright triples. We conclude 
that optically-faint triples would not populate the lower- 
right corner of panel A, and that the absence of such a 
population represents the true properties of quasars. 

The absence of quasars in the upper-left corner of 
panel A (with high R and low Ri) is also due to in- 
trinsic quasar properties, as demonstrated in panel C. 
Quasars can be classified as triples only if they have two 
lobes with S'lobo > ImJy; otherwise, one or both lobes 
would be undetected. Because of the optical and radio 
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Fig. 5. — Number of spectra contributing to each spectral com- 
posite as a function of wavelength bin. Colors correspond to the 
colors of spectral composites in Figure [G] 

detection limits of our sample, quasars with high R and 
low if they exist, are likely to have (an) undetected 
lobe(s), and therefore be classified as core or lobe rather 
than triple. Panel C compares the Ri distribution of the 
triples to that of the lobe plus core classes. The latter 
do not show a lower range of which indicates that 
no such population of quasars with high R and low Ri 
exists. 

In summary, Figure H] shows a strong correlation be- 
tween R and We demonstrated, using distributions 
in Panels B and C, that the observed correlation is 
not a selection effect, but is instead an intrinsic prop- 
erty: quasars with high R tend to have high and 
vice-versa. This observation supports the hypothesis 
that both parameters are measures of quasar orienta- 
tion. The scatter in the correlation supports the idea 
that other factors are also influencing these two mea- 
surements. 

5. THE DEPENDENCE OF OPTICAL SPECTRA ON 
RADIO SELECTION 

The spectral composite is a valuable tool for 
studying properties of large samples of astronomi- 
cal sources. In recent years, the advent of large 
spectroscopic surveys has allowed for in creasingly de- 
tailed s tudies of AGN spect r al com p osites (iFrancis et al.l 
1991; 'Baker & Hunstcad TogF; 'Zheng et alJ |1997|; 
Brothcrton ct al. 2001; Vanden Berk ct al. 20Ql5. Spec- 
tra of many objects can be combined to create compos- 
ites with high signal-to-noise, allowing average prop- 
erties to be compared across sample types. Here we 
compare spectral properties of quasar radio morphol- 
ogy classes, and investigate whether a selection based 
on radio properties implies concomitant changes in op- 
tical spectra. 

5.1. Construction of quasar spectral composites 

Creating a spectral composite involves several steps: 
shifting each input spectrum to its rest frame, rebinning 



the spectra to matching wavelength grids and resolution, 
normalizing their fluxes, and stacking into a final com- 
posite. Variations of these steps can lead to significant 
differ ences in the resulting composite ([Francis et al.l 

[TMI). 

We redshift-corrected each spectrum using the S07 
redshift fSection l2.2p . Spectra were then re-binned, con- 
serving flux, onto a common wavelength grid with bins 
~ 69 km s~^ wide, the same sampling as observed SDSS 
spectra. 

We used the iVanden Berk et al.l ()2001[ ) quasar spec- 
tral composite as a template for rescaling our input spec- 
tra to a common flux level, using the following method. 
In the region of wavelength overlap, we found the flux 
ratio of the input spectrum to the Vanden Berk et al. 
composite at each wavelength; however we ignored the 
region blueward of the Lya emission line (A < 1250A) so 
as to avoid uncertainties due to redshift-dependent ab- 
sorption in the Lya forest region. We then determined 
the median flux ratio in the overlap region, and used 
this value to rescale the input spectrum. Our resulting 
composites are in arbitrary flux density units (but on 
the same scale as the Vanden Berk et al. composite). 

We stacked the rescaled spectra into two different fi- 
nal composites: a "median" composite and a "geometric 
mean" composite. (For the latter, negative flux val- 
ues were ignored.) A median composite preserves the 
relative fluxes of emission features, while a geometric 
mean preser ves the overall continuum shape for power- 
law spectra (jVanden Berk et "all 1200 1[ ). We note that 
our composites do not include BAL quasars, as BALs 
can influence a spectrum's shape. 

We compared the median composite for the 72,223 
quasars (both radio and radio-quiet) without BALs to 
the median composite from Vanden Berk et al., which 
also excludes BAL quasars. The two spectra agree very 
well and are nearly indistinguishable to the eye, except 
at the extreme wavelength values, where fewer spectra 
contribute and thus the composites are noisier. 

The spectral composites for the entire sample and for 
individual morphological classes are discussed in the re- 
mainder of this section. Table [5] presents the composite 
for the radio-quiet quasars, for the full sample of radio 
quasars, and for the radio morphology classes. 

Individual quasar spectra have different rest-frame 
wavelength coverage owing to their wide redshift range 
(0.78 < z < 5.4). For each spectral composite, differ- 
ent objects contribute to the long-wavelength and short- 
wavelength regions. Figure [5] shows the number of input 
spectra that make up each section of a spectral compos- 
ite as a function of wavelength. 

5.2. Spectral composites of individual radio morphology 

classes 

Geometric mean composites for distinct radio mor- 
phology classes are compared by ratio in Figure [HI Dif- 
ferences in individual composites are difficult to distin- 
guish by eye, but are emphasized in ratios of the compos- 
ites. We have separated the core sources into resolved 
and unresolved ( §3.2.3p . Because those two groups dom- 
inate the radio quasars, the core/radio ratios are nearly 
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TABLE 5 

Median Composite Quasar Spectra for Radio Morphology Classes 



A 


Radio-quiet 


All Radio 


Res. 


Core 


Unres. 


Core 


Jet 




Lobe 


Triple 


(A) 


fx 




fx 




fx 




fx 




fx 




fx 




fx 




900.00 


1.80 


0.148 


1.90 


0.681 


0.664 


1.65 


2.18 


0.822 


1.39 


3.78 


2.66 


1.34 


4.40 


-1.00 


900.21 


1.86 


0.149 


2.19 


0.693 


0.980 


1.58 


2.47 


0.844 


1.83 


1.80 


2.22 


1.78 


5.37 


-1.00 


900.41 


1.77 


0.144 


1.80 


0.713 


0.808 


1.38 


1.81 


0.847 


0.488 


4.16 


3.57 


2.68 


4.96 


-1.00 


900.62 


1.90 


0.143 


1.76 


0.740 


0.880 


2.40 


2.09 


0.847 


1.13 


3.01 


1.19 


2.83 


2.28 


-1.00 


900.83 


1.88 


0.140 


1.95 


0.798 


1.28 


1.75 


1.78 


0.921 


2.83 


2.76 


2.06 


4.22 


3.79 


-1.00 



Note. — The full table is available in the electro nic version of the paper. Values of fx (f lux) are in arbitrary 
units, such that they are on the same flux scale as the lVanden Berk et al.l 1120011) spectrum (see ^S.ll for details). The 
uncertainty in the median is given hj af, a value of —1 indicates that only one input spectrum contributes to that 
wavelength bin, while indicates that no input spectra contribute at that wavelength. 



Median Composite 



TABLE 6 

)UASAR Spectra for Triple 



JuASAR LoG(i?) Subclasses 



A 

(A) 


log(R) < 
fx 


-1.5 


-1.5 < log(i?,) < -0.5 
fx cr/ 


-0.5 
fx 


< log(ij) < 0.5 


0.5 < log(R) 
fx o-/ 


920.14 


24.3 


1.00 


0.00 


0.00 


0.00 


0.00 


0.00 0.00 


920.35 


19.3 


1.00 


0.00 


0.00 


0.00 


0.00 


0.00 0.00 


920.56 


6.94 


1.00 


0.00 


0.00 


0.00 


0.00 


0.00 0.00 


920.78 


8.78 


1.00 


0.00 


0.00 


0.00 


0.00 


0.00 0.00 


920.99 


6.27 


1.00 


0.00 


0.00 


0.00 


0.00 


0.00 0.00 



Note. 



The full table is available in the electronic version of this paper. Values of 



are in arbitrary units, such that they are on the same flux scale as the lVanden Berk et al.l 11200 II ) 
spectrum (see ^5.11 for details). The uncertainty in the median is given by tr^; a value of —1 
indicates that only one input spectrum contributes to that wavelength bin, while indicates that 
no input spectra contribute at that wavelength. 



equal to 1 at all wavelengths. Composites for different 
morphology classes show similar spectral shapes. The 
second panel suggests that the average optical spectral 
index is redder for radio quasars than for radio-quiet 
quasars (see §3.4l for further discussion). However, spec- 
tral indices measured from composites may not be accu- 
rate because objects at different redshifts contribute at 
different wavelengths owing to varied rest-frame wave- 
length coverage. Thus the spectral index of a composite 
depends on quasar evolution. 

The ratios of spectral composites show that the 
strength of spectral lines vary with radio power and 
morphology. For example, CIV A1549 is strongest in 
jet and triple sources. Figure [7] directly compares the 
line profiles of the median spectral composites. Spectra 
in this figure have been scaled such that the continuum 
value at the line center is equal to 1. 

In Figure El the [OIII] doublet and [Oil] are clearly 
stronger in radio quasars than radio-quiet quasars, an 
effect that is expected from previous o bservations of 
these lines. iRawlings fc SaundersI ()1991[ ) observed that 
the total narrow line luminosity in quasars (measured 
from [Oil] and [OIII]) positively correlates with jet 
power, as measured from the lobe emission and spec- 
tral ages. In our sample, these lines are strongest in the 
triples and weakest in the lobes and jets. These results 
agree with Rawlings & Saunders, given that triples have 
greater extended emission than the lobe and jet sources 
(which suggests that triples have the most powerful ra- 



dio jets). The resolved cores are weaker in [OIII] than 
the triples, but stronger than the lobes and jets, and are 
even stronger in [Oil] than the triples. This evidence 
suggests that the resolved cores contain many powerful 
sources. As surmised in §3.21 these sources are probably 
still young, and therefore small in physical extent, which 
explains why their lobes are not resolved into separate 
components in FIRST. 

Among the broad lines, CIII] and Mgll show no sig- 
nificant change among the radio classes. However, Lya 
and C IV show strong variation, and are clearly the 
strongest in the triple and jet samples. 

5.3. The effects of orientation on composite spectra 

By dividing the sample of triples according to R, we 
can create spectral composites for quasars with different 
ranges of assumed o rientation angle to the radio jet axis. 
iBaker fc HunsteadI ([l995, ) performed a similar analysis 
of 47 radio quasars. Their results suggest that as the 
viewing angle to the jet axis increases, the optical con- 
tinuum steepens and line widths increase. We now re- 
visit these findings using an order-of-magnitude larger 
sample. 

We divided the triples into the following four groups 
based on the rest-frame core-to-lobe ratio R at 20 cm 
(Eq. [31): (1) 0.5 < log(i?), (2) -0.5 < log(i?) < 0.5, 
(3) -1.5 < log(i?) < -0.5, (4) log(i?) < -1.5. These 
groups contain 30, 277, 190, and 20 triples, respectively. 
Groups 2, 3, and 4 roughly correspond to the three 
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Fig. 6. — Comparison of geometric mean spectral composites 
of radio quasars according to morphology. Top panel: composite 
spectrum for radio quasars. Second panel; ratio of radio quasar 
composite to composite for radio-quiet quasars. Third panel; ratio 
of composite for unresolved core sources to all radio sources. The 
remaining four panels show analogous ratios for resolved cores, 
jets, lobes, and triples compared to the full radio sample. 



groups defined bv lBaker &: HunsteadI ()1995[ ). who used a 
measure of i? at 3 cm in the rest frame. The conversion 
between the 3 cm and 20 cm values is 



R[20cm] =(20/3)-W5) xR[3„„]. 



(8) 



Based on the quasar core-boosting models, group 1 
should correspond to quasars with lines-of-sight close 
to the radio-jet axis, while group 4 should contain 
quasars viewed closest to the plane of the accretion disk. 
Group 4 contains only 20 sources from the main radio 
quasar sample. We add to this category by including 
double-lobed quasars, which are clearly lobe-dominated 
by virtue of having no radio-detected core. We calcu- 
lated an upper limit on R for the double-lobed quasars, 
assuming a core flux density less than the FIRST sur- 
vey limit of 1 mjy. Ninety-seven double-lobed quasars 
have an upper limit of log(i?) < —1.5, and we included 
their spectra when making the composite spectrum for 
group 4. The median spectral composites for the four 
log(i?) groups are presented in Table [H 
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Fig. 7. — Spectral line profiles from the median composite spec- 
tra of each radio morphology class. Colors correspond to the colors 
in Figure [5] cyan — unresolved core; purple — resolved core; blue, 
jet; green, lobe; red, triple. The spectral line profiles were fit to a 
simple Gaussian -|- line continuum model, and have been rescaled 
such that the continuum at the line center is equal to 1. 



The geometric mean composites are shown in Fig- 
ure [51 We compare them by examining the ratio of 
each with the compo site spectrum for al l triple s. Con- 
trary to the results of IBaker &: HunsteadI (jl995D . we see 
no change in spectral shape with R. However, there is 
a clear change in the CIV A1549 line profile: EWciv 
appears to increase monotonically toward sources with 
the highest R. 

Individual line profiles are shown in Figure |9] us- 
ing the median composites. It appears that the 
EW of the broad lines (CIV, CIII], Mgll and the 
Balmer lines) increase with R. In contrast, the 
[OH] and [OIII] lines show roughly the opposite 
trend: their EW appears to be anti-correlated with 
R. Previous observations have suggested that both 
narrow and broad line E W are anti- correlate d wit h 
R, as claimed for Mgll b v lBrowne h Murphvl ()1987| ). 
IBaker fc Hunste ad' (1995'), and 'Baker (1997); for the 
Balmer li nes by Baker & Hunstea d (1995.) : and for 



a d 

O HI] by iWills fc Brownei (|1986D , ljacksOTL' fc Brow^ 
IMl) . and lBaked (11991 7" These results from our com- 



posite spectra are consistent with previous claims for 
the narrow lines, but are inconsistent for the broad lines. 
However, we note that the scatter of EW in individual 
spectra contributing to each composite is large, such 
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Fig. 8. — Comparison of geometric mean spectral composites 
for quasars categorized by core-to-lobe ratio at 20 cm in the rest 
frame. Top panel: composite spectrum for all triples. Second 
panel from top; ratio of triple composite to radio quasar compos- 
ite. Remaining panels: ratios of other spectral composites com- 
pared to the spectrum shown in the top panel. The third, fourth, 
and fifth panels (blue, red, and green) correspond to subsets of the 
triples sample. The bottom panel includes double-lobed quasars 
without a radio core. 



that the uncertainty in the median EW for each hne 
tends to be much greater than the change in EW seen 
from one composite spectrum to the next. We will in- 
vestigate these trends quantitatively in fj6l the results of 
correlation analyses presented there suggest that emis- 
sion line EW is not significantly correlated with R. 

Motivated by our conclusion in Section U that R and 
i?i are physically correlated, and by evidence that the 
latter is the stronger statistical measure of o rientation 
(| Wills fc BrothertonI [19951 IKharb et aII[20lol) . we now 
investigate composite spectra grouped by Ri. We cre- 
ated the following four groups: (1) 2.5 < log(i?), (2) 
2.0 < log(i?) < 2.5, (3) 1.5 < log(i?) < 2.0, (4) log(i?) 
< 1.5. These groups contain 1070, 960, 1412, and 1272 
quasars, respectively. These groups are much larger 
than the log(i?) groups because log(i?i) can be deter- 
mined for every quasar, not just for triples. 

The geometric mean composites for the four Ri quasar 
groups are shown in Figure fTUl Strong differences in Lya 
and C IV are apparent in the ratios of composites, as 
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Fig. 9. — Spectral line profiles from the median composite spec- 
tra of quasars categorized by core-to-lobc ratio at 20 cm in the 
rest frame. The samples and colors in this figure correspond to 
those in FigurefS] black, all triples; blue, triples with log(/?) > 0.5; 
red, triples with —0.5 < log(ij) < 0.5; green, triples with —1.5 < 
log(i?) < —0.5; orange, triples with log(_R) < —1.5. The spectral 
line profiles were fit to a simple Gaussian -|- line continuum model, 
and have been rescaled such that the continuum at the line center 
is equal to 1. 

well as subtle differences in CIII], Mgll, and H/3. Each 
of these lines appears to be the strongest in group (1), 
with decreasing strength toward group (4). These differ- 
ences are even more clear in Figure lll[ which shows the 
emission lines in the median spectral composites. This 
figure indicates that the EW of quasar emission lines 
increases as i?i increases. In the following section, we 
investigate these intriguing trends more quantitatively 
by examining correlations among individual spectra. 1 

6. CORRELATIONS OF SPECTRAL LINE PARAMETERS 
WITH R AND Ri 

Three commonly- accepted theories that relate quasar 
line behavior to viewing angle are: orientation- 
dependent obscurat ion by a clumpy dust torus (e.g., 
iNenkova et al. 2008), Dopplcr-boosting of the contin- 
uum (e.g., Browne & Murp hy .1987) , a nd inc lination 
angle of the accretion disk iNetzeii le.g- fl987l) . Each 
model predicts that the EW of quasar emission lines 
increases as the viewing angle to the radio-jet axis in- 
creases. These theories share the assumptions that the 
line emission is isotropic while the continuum emission 
depends on the orientation to the accretion disk. 
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TABLE 7 

Median Composite Quasar Spectra for. Quasar Log(/?i) Subclasses 



A 


log(fii) < 1.5 


1.5 < log(i?.i) < 2.0 


2.0 < log(i?,i) < 2.5 


2.5 < log(i?i) 


(A) 


fx 




fx 




fx 




fx o-/ 


900.00 


3.18 


0.814 


1.79 


1.36 


1.48 


1.27 


1.90 1.76 


900.21 


1.88 


0.938 


2.19 


1.12 


1.83 


1.48 


2.90 1.93 


900.41 


0.956 


1.17 


1.67 


1.36 


1.63 


1.05 


2.07 2.05 


900.62 


1.25 


1.11 


2.83 


1.30 


2.02 


1.50 


1.28 1.96 


900.83 


1.29 


1.24 


2.79 


1.42 


2.65 


1.18 


1.25 2.55 



Note. — The full table is available in the electronic version of this pap er. Values of fx (flux) 
are in arbitrary unit s, such that they are on the same flux scale as the Va nden Berk et al.l BOOl. ) 
spectrum (see i|5.1l for details). The uncertainty in the median is given by trj; a value of — 1 
indicates that only one input spectrum contributes to that wavelength bin, while indicates that 
no input spectra contribute at that wavelength. 
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Fig. 10. — Comparison of geometric mean spectral composites 
for quasars categorized by log(iJi). The top panel shows the com- 
posite spectrum for all radio quasars. The remaining panels each 
show the ratio of a composite for a subset of quasars (as labeled) 
to the composite in the top panel. 

In this section, we examine the dependence of quasar 
emission lines on the core-boosting parameters R and 
We find no significant correlations with i?, but find 
several significant correlations with Ri in the opposite 
direction from what is predicted by the above models. 

6.1. Measuring correlations in the data 

To evaluate correlations between line par amete r and R 
or we ran the linear regression code of lKellvl (|2007l ). 
which uses Bayesian methods to incorporate uncertain- 
ties in the data, as well as upper-limits. The parame- 
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Fig. 11. — Spectral line profiles from the median composite spec- 
tra of radio quasars categorized by core radio-to-optical flux den- 
sity ratio in the rest frame. The samples and colors correspond 
to those in Figure [TUl blue, log(iJi)> 2.5; red, 2 < log(i?i) < 2.5; 
green, 1.5 < log(iJi) < 2; orange, log(i?i) < 1.5. The spectral line 
profiles were fit to a simple Gaussian + line continuum model, 
and have been rescaled such that the continuum at the line center 
is equal to 1. 



terization we used is 

log(r) 



A + Bx log(X), 



(9) 



where Y is the dependent variable (e.g., EW), X is the 
independent variable (e.g., R or and A and B are 
intercept and slope respectively. All line parameters 
X have 1-cr error bars ax determined from lShen et alj 
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((Moh . For sources where X — ax < 0, we treated the 
source as a non-detection with upper- Hmit oi X + 3crx • 
This modification was necessary in approximately 1.5% 
of sources. 

To investigate correlations with R, we used the sample 
of 519 triple quasars and 383 lobe quasars with known 
core and lobe flux densities (i j3.2.4l)F^ To investigate 
correlations with we used the full radio quasar sam- 
ple. We limited the analysis to spectra with median 
signal-to-noise per pixel (SDSS pixels are 69 km s~^ 
wide) of at least 5 in the relevant line region (|Shen et al.l 
IMm . We excluded sources with a bro a.d absorption 
line fe ature near the line, as identified by iGibson et al.l 
(|2009D ; such features can affect the measurement of line 
EW. We also excluded sources with the line in absorp- 
tion. Such exclusions do not hinde r the B ayesian linear 
regression analysis, as discussed in iKellvl (2007). 

In addition to estimating the intercept A and slope i?, 
the code estimates the Pearson correlation coefficient. 
The latter ranges from —1 to 1 where indicates no 
correlation, 1 indicates a perfect correlation, and — 1 
indicates a perfect anti-correlation. We consider a cor- 
relation to be significant if both the slope B and the 
correlation coefficient are non-zero with certainty > 3(T. 
Results of each test are listed in Table [SJ and are dis- 
cussed qualitatively in the remainder of this section. 

The scatter in EW measurements is much larger than 
errors in individual data points. A number of intrin- 
sic processes may contribute to the scatter. Quasars 
are known to vary in the optical co ntinuum by abou t 
20% over multi-year timescales (e.g.. lHook et al.iri994[ ). 
which may be due to changes in the accretion disk or jet 
variation s (iMarscher &: Geaii 119851) . As demonstrated 
by, e.g., iVanden Berk et aTT (|2004[ ) . variability ampli- 
tude in the optical depends on wavelength (quasars 
are more variable at shorter optical wavelengths), lu- 
minosity (more luminous quasars vary more strongly), 
and redshift (higher redshift quasars show stronger 
variability). Radio variab ility of the core can occur 
at the 20% level (Barvain is et al.l l2005| ) ahhough this 
is typically observed in only 1-2% of FIRS T sources 
(|R,vs fc Machalskil \Tmi Ide Vries et all [200l . Differ- 
ences in local environment may also cause varia tions in 
emission line luminosities (e.g.. lMcCarthviri993[ ). 

6.2. Equivalent width of Mg II and C IV broad emission 

lines 

We begin by investigating the EW of Mgll A2799 
(EWMgii) as a function of R and Ri (Fig. The 
left panel of Figure [T^] shows EWivig ii vs R for triple- 
morphology quasars, using the average lobe flux density 
to calculate R. The middle panel shows EWMgii vs R 
for /o&e-morphology quasars, using the single detected 
lobe to calculate R. For plotting purposes only, we show 
median values of EWMgii in Alog(i?) — 0.5 bins to help 
guide the eye. To perform the linear regression analysis, 
we used the true data points, not the median values. 

The remaining 100 triple quasars and 4 lobe quasars were 
excluded due to crowded FIRST images. The identities of their 
l obe co mponents, and therefore their R values, are ambiguous 

(E23J- 



There is no significant correlation of EWivig ii with i?; 
this result is reflected quantitatively via the results of 
the linear regression (Table [H). However, we do see a 
strong correlation with Ri in the rightmost panel. In 
the full sample of radio quasars, the correlation is sig- 
nificant at approximately the lOtr level. When divid- 
ing the radio quasars sample into individual morphol- 
ogy classes, we continue to see a significant correlation, 
but with reduced strength (~4-7ct) in the core, triple, 
and jet classes. The results of the linear regression 
vary across the different morphology classes; the triple 
and jet classes have the steepest slopes while the core- 
morphology sources show weaker, but still significant, 
correlations. The lobe sources show the same trend, 
but not at a significant level (only l.Gcr). The greatest 
variation is seen between the unresolved core class and 
the triple class, which differ in slope by 3cr. We con- 
clude that there is certainly a positive trend of EWivigii 
with Ri in radio quasars, and that the trend may vary 
weakly with morphology class. 

Figure [13] shows the equivalent analysis for EWc iv • 
Similar to the results for Mgll, we see no correlation 
of EW with i?, but we do see a significant positive cor- 
relation (~ Sct) with Ri. The same trend with Ri is 
observed in the individual morphology classes, but at 
a weaker level. We observe significant correlations for 
the core and triple classes, but not for the lobe or jet 
classes. In this case, the overall variation among the 
different classes is only ^ 1.2a. 

Our observations conflict with the results of iBaked 
()1997D . Baker observed a significant anti-correlation of 
EWMgii with i? in a sample of 67 quasars with z < 2.2; 
a similar qualitative, but not significant, trend was ob- 
served for EWc IV- Baker reported a stronger trend for 
EWMgii in quasars with z > 1 than in those with z < 1. 
We see no difference in correlations with R when testing 
samples in these two redshift ranges. Baker proposed 
orientation-dependent obscuration by a dusty torus as 
an explanation for their observed trend, suggesting that 
the torus increasingly obscures the accretion disk at high 
inclination angles. We note that, while our quasars 
match the Baker sources in redshift (e.g., z < 2.2 for 
Mgll emitters), our samples differ in radio brightness: 
the Baker sources had radio flux densities S > 0.95 Jy 
at 408 MHz; the brightest Mgll triple in our sample 
has a radio core flux density of 0.71 Jy. If orientation- 
dependent obscuration is responsible for the trend ob- 
served by Baker, our results may show that dust obscu- 
ration does not depend on orientation in these lower- 
luminosity quasars. Alternatively, our results could in- 
dicate that the continuum emission and the broad line 
region (BLR) intrinsically have the same dependence on 
orientation such that EW does not change with orienta- 
tion. This situation could occur if both the broad line 
emission and the c ontinuum emission o riginate in the 
accretion disk (e.g.. lCollin-Souffrinlll987D . 

Why is there a significant correlation of EW with Ri 
but not with R.7 The Ri distribution may have a higher 
signal-to-noise ratio because it includes a larger number 
of data points: R can only be determined for quasars 
with at least one observed radio lobe, while Ri can be 
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Fig. 12. — EW of Mgll A2799 as a function of core-boosting parameters R and ijj. In the left and middle panels, the full distributions 
are shown as black points. In the right panel, the full distribution is shown as linearly-spaced contours. The number of sources in each 
panel and each morphology class is listed in Table |8] To guide the eye, colored error bars show the median (and error in the median) EW 
in evenly-spaced log(i?) or log(iJi) bins. Horizontal dotted lines show the median value for each panel. Dashed lines show the best-fit 
line for the full set of sources in each panel. Quantitative results of the linear regression for all morphology classes are listed in Table |8] 
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Fig. 13. — EW of CIV A1549 as a function of core-boosting parameters R and Ri. Colors and symbols are as in Figure [T2l Number of 
sources of each class and results of the linear regression analysis are listed in Table [8] 



determined for all radio quasars. However, we observe 
a stronger Ri dependence even among the triple sources 
alone. One possible explanation is that the EW-i?i 
correlation is based on some intrinsic physical cause, 
rather than on orientation. However, if both param- 
eters are statistical indicators of orientation, then our 
data suggest that broad line EW does depend on ori- 
entation, and that i?i is the stronger statistical orienta- 
tion param eter. This conclusion is co nsist ent with the 
propo sal of lWills fc BrothertonI (|1995D and lKharb et al.l 
(|2010D that R additionally depends on local environ- 
ment, while i?i does not. One possible explanation for 
these i?i-EW correlations is an anisotropically-emitting 
BLR, composed of spherical clouds that emit most in- 
tensely from the inner side of each cloud, illuminated 



by the ionizing dis k emission fe.g.. iFerland et aDll992l : 
'Koris ta et al.l[l997l ) . An edge-on observer would see ion- 
ized emission only from clou ds on the far side of the BLR 
for, e.g., CIV (see Fig. 4 of iHorne et aL|[200l . In con- 
trast, a face-on observer would see half the ionized face 
of every BLR cloud, resulting in more BLR emission 
relative to the continuum and therefore a larger EW. 

6.3. Equivalent width of the [O III] narrow line 

FigurelTHshows the distribution of EW[oiii] as a func- 
tion of R and i?i . We do not observe any significant cor- 
relations with R nor with i?i . We see a trend in the full 
radio sample for EW[oiii] to increase with but with 
a significance of only l.Scr; no individual morphology 
classes show any stronger trend. 
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Fig. 14. — EW of [OIII] AA4959/5007 as a function of core-boosting parameters R and Ri. Colors and symbols are as in Figure [T2l 
Number of sources of each class and results of the linear regression analysis are listed in Table \E\ 

The Baldwin effect is a well-known anti-correlation 
between broad-line width and optical continuum lum i- 
nosity. It was first observed in CIV (jBaldwinl I1977D , 
and has since been disco vered in many other lines (e.g., 
lOsmer fc Shields 1999; Wu et al.l 120091 and references 
therein) . In this section, we compare the Baldwin effect 
in our sample to previous observations. 

Traditionally, the Baldwin effect has been measured as 
EWciv vs monochromatic luminosity at 1450 A. Here 
we use the continuum level at the rest-wavelength of the 
CIV line, 1549AF1. This modification should not signif- 
icantly affect the slope of our Baldwin effect measure- 
ment because the lumin osities at these two wavelength s 
are strongly correlated (|Yee fc Okelll978HShudellll98lD . 

Figure [15] shows the Baldwin effect we observe in each 
quasar morphology class. Previously observed slopes 
are t ypically around —0.23 (jOsmer et al.l 1199^ IGreenI 
119961 e.g.). Our radio-quiet sample has a slightly steeper 
slope of B = -0.257 ± 0.005. The slopes for the radio 
quasars are consistent with this value, ranging from B = 
-0.180 ± 0.085 to B = -0.297 ± 0.068 for ah classes 
except jet sources. The jet class, with a slope of i? = 
—0.098 ± 0.114, has a slope consistent with zero. 

It has been noted that the Baldwin effect appears to 
be shallower for lower-lu minosity quasars. For exam- 
ple, [Dietrich et al.l ()2002f ) observed a shallower slope of 
—0.14 ± 0.02 in a sample of quasars spanning six or- 
ders of magnitude in luminosity, but observed a slope 
of —0.20 ± 0.03 when limiting their sample to quasars 
with_Ml450A) > 5 X lO^^ erg s^^ Hz^i (but see also 
iWu et al.|[2009l ). For most of our sources, our slope val- 
ues are consis t ent w ith the higher-luminosity slope of 
iDietrich et aD (|200l : all of our quasars are brighter 
than the luminosity limit they used to separate their 
sample into high- and low-luminosity samples. 




29.5 



30.0 30.5 31.0 

log(L,54<,A[erg/s/Hz]) 



31.5 



Fig. 15. — The Baldwin effect across the quasar radio morphol- 
ogy classes: EWc iv as a function of continuum luminosity at 1549 
A. Symbols correspond to median values in bins with a width of 
0.4 in log(L^) from 29.5 to 31.5. 



The weak positive correlation between narrow line 
EW an d core-boosti ng p arameter conflicts wi th the re- 
sults of lBakeii (|1997D and lJackson et al.l ()1989[ ). who ob- 
served an anti-correlation of EW[oiii] with R. We note 
that the quasar samples used in those studies were an 
order of magnitude smaller than the sample presented 
in this paper. As for Mgll reported above, our sample 
matches the Baker sample in redshift, but not in radio 
brightness (our sources are much fainter). While the 
orientation-based theories discussed in the beginning of 
Section [5] suggest that EW[oiii] is anti-correlated with 
core-boosting, we see no evidence of such a trend in our 
data. 



7. SUMMARY AND CONCLUSIONS 



6.4. The CIV XI549 Baldwin Effect 



Continuum measurements at 1549A are from the SDSS. 



18 



Kimball et al. 



We have compiled the largest sample to date of ra- 
dio quasars with high-quality spectra and reliable visual 
classification by radio morphology. The sample consists 
of 4714 radio quasars from the SDSS DR5 spectroscopic 
survey with a 20 cm FIRST detection brighter than 
2 mJy within 2" of the optical position. We have used 
the sample to investigate how optical spectral proper- 
ties of a quasar sample depend on a purely radio-based 
selection. One investigative technique involves spec- 
tral composites generated for individual radio morphol- 
ogy classes. We also compared two measures of quasar 
core-boosting, R and and searched for correlations 
between these parameters and spectral line equivalent 
widths. Data products available in the electronic ver- 
sion of the paper are listed in the Appendix. 

Our main results are as follows: 

1. We found that generally it is not true that radio 
quasars have redder optical colors than radio-quiet 
quasars. Among a z < 2.7 sample, radio quasars 
with unresolved core morphology contain a higher 
fraction of extremely-reddened objects than do ra- 
dio quasars with resolved core, triple, lobe, jet, 
or double-lobe morphology, as well as radio-quiet 
quasars. These results also suggest that the SDSS 
color-selection criteria for targeting quasars are 
not strongly (< 0.1 mag) biased, because all other 
radio quasar classes have observed color distribu- 
tions similar to that of radio-quiet quasars. 

2. We found that Ri and R, two parameters that 
have been proposed as measures of quasar radio- 
core boosting, are correlated. This observation in- 
dicates that some physical factor drives both pa- 
rameters; core-boosting models suggests that ori- 
entation is the driving factor. The scatter in the 
observed correlation indicates that another effect 
is also involved; external environment or age are 
two such possible factors. 

3. In contrast to iBakeii ()1997D . we did not find 
significant correlations between R and emission 
line EW for CIV A1549, Mgll A2799, or [OIII] 
AA4959/5007. These results suggest that orienta- 
tion does not govern EW of quasar emission lines, 
or that the line emission and the optical contin- 
uum have the same dependence on orientation, 
as might be the case if both the lines and con- 
tinuum originate in a disk-like geometry. On the 
other hand, we did observe positive correlations 
between i?i and broad line EW for the full radio 
quasar sample and also for some of the individ- 
ual morphology classes, including core- and triple- 
morphology sources. If Ri is taken as a statistical 
indicator of orientation, then these observations 
suggest that orientation does play a role in deter- 
mining broad line EW in radio quasars. 

At first glance, the above conclusions appear to 
be inconsistent, given that both R and Ri should 
corre late with line-of-sight or ienta tion. How- 
ever, iWills fc Brothertoij (jl995[ ) and IKharb et ahl 



(|2010[ ) suggested that i?i is a stronger orienta- 
tion indicator than R, because the latter parame- 
ter also depends significantly on environment and 
age, as discussed in Section S] Our results are 
consistent with the idea that line EW depends on 
orientation, because of the observed correlations 
with Ri. A corollary of that conclusion is a pos- 
sible confirmation that Ri depends more strongly 
on orientation than does R, because of the lack 
of observed correlations between EW and R. In 
that case, the observed correlations may be the 
result of an anisotropically-emitting BLR (e.g., 
iHorne et al.l [2003i) . However, our results do not 
constitute definitive proof of these conclusions. It 
is also possible that a physical parameter other 
than orientation is the driving factor of the Ri- 
EW correlations. 

4. We observed the Baldwin effect (the anti- 
correlation between EWciv and continuum lu- 
minosity) in both radio and radio-quiet quasars. 
We find, for most morphology classes, slopes of 
approximately —0.25. Owing to large uncertain- 
ties, the slope for the double-lobed classes is not 
strong enough to be considered a significant trend, 
while the slope for the jet classes is consistent with 
zero. For all other morphology classes (the vast 
majority of the sample), the observed trends are 
significant; these results are generally consistent 
with previous observations of the Baldwin effect 
in high-luminosity quasar samples. 

If Ri is a statistical indicator of orientation as pro- 
posed, our results suggest a change in quasar emis- 
sion line EW with viewing angle. For isotropic 
line emission, the observed correlations imply that 
the continuum emission increases as the viewing an- 
gle to the radio-jet axis increases. This result is 
the opposite of what has been o bserved in earlier, 
much s maller, s amples of quasars: 'Browne fc Mu rphvl 
(fT987l) ■ [BakeFfc Hunstead (19951 and Baker (1991) ob- 
served an anti- correlation between c ore-d ominance and 
EWMgii, while I Jackson et al.l (|1989D and|Baker| ()1997f ) 
observe a similar anti-correlation among the narrow 
emission lines. Basic core-boosting and disk-inclination 
models predict an increase in continuum emission for 
small viewing angles (with respect to the radio-jet 
axis), which would lead to smaller equivalent widths 
for isotropically emitted lines from quasars with high 
core-dominance. The evidence for increased broad line 
EW with Ri in this sample is strong. It is possible 
that the interpretation of Ri as an orientation indicator 
is not correct, although its correlation with R argues 
against this conclusion. Our results may indicate, there- 
fore, that the broad lines are not emitted isotropically. 
A significantly flattened broad line region (BLR) might 
produce the observed trends with Ri. In that case, these 
trends would imply a quite anisotropic distribution of 
BLR clouds, or perhaps a BLR that is associated with 
the accretion disk. 

We extend thanks to Patrick Hall, Brandon Kelly, Rob 
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Linear regression results in Section [6] 
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Note. — Distributions with significant correlations are shown in bold text. We define a significant correlation as 
having slope B and correlation coeflicient that are non-zero with at least 3(t certainty. We exclude statistical analysis 
for samples with fewer than 20 sources. 

t Parameters given in Equation [9] 

t Value of Pearson correlation coefficient reported by the code of IKellvl H2007I V 
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APPENDIX 

OPTICALLY-FAINT TRIPLE-MORPHOLOGY QUASARS 

We discuss here the construction of the sample of "optically-faint triple-morphology quasars" : triple-morphology 
sources with no optical detection, or with an optical source too faint for the SDSS spectroscopic quasar sample. We 
compile this sample in order to determine whether optically-faint triples have a similar R distribution to optically- 
bright triples. The comparison of the two distributions represents one step toward demonstrating that the observed 
correlation between R and Ri is an intrinsic property of radio quasars ([j4|). 

We drew candidat es from the ra dio catalog of KI08, which identifies all sources from the 20 cm NRAO-VLA 
Sky Survey (NVSS) HCondon et al.l lT998) with three FIRST matches within 120". (Only one triple quasar in our 
main sample is undetected by NVSS.) The initial selection criteria, based on parameters in the KI08 catalog, are: 
matchf lag_nvss = —1, matchf lagjirst = 1, matchf lag_120 = 3. These criteria select unique NVSS sources with 
3 FIRST matches within 120". 

Using simple geometry, we eliminated all configurations where the FIRST triplet is arranged in an acute triangle 
on the sky. Such configurations are almost certainly spurious: extremely few double-lobed quasars have an opening 
angle as small as 90°, as observed by dV06 and confirmed in our main sample triples fFigure fT6|) . We identify the 
core component as the one "in-between" the other two FIRST sources. Following our main sample selection, we 
require 5*0010 > 2 mJy. To further remove spurious sources, we require that the core source be unresolved in the 
FIRST survey: S'20/5'pcak < 1-23 (KI08). 

We then eliminated configurations with an optical core bright enough for SDSS targeting: we matched the cores 
to DR6 within 2", and removed any with an optical match brighter than i — 19.1, as well as others which already 
had an SDSS spectrum. 

We decided the final selection criteria by comparing parameters of the remaining 4361 candidates to the main 
sample of triples. Figure [16] shows the lobe — core — lobe opening angle and lobe — lobe magnitude differences (fiux 
density ratios) for the "optically-faint" and "optically-bright" triples. Based on the distribution of the optically- 
bright sample, we require that the optically-faint candidates have an opening angle of at least 120° and a lobe — lobe 
flux density ratio between 0.1 and 10 (corresponding to a magnitude difference of less than 2.5 mag). The final 
sample of optically-faint triples contains 2700 sources. 

Of the 2700 optically-faint triples, 719 have an SDSS object located within 2" of the core. These sources, listed in 
Table [9l are excellent candidates for optical quasars, but were too faint for SDSS spectroscopy. 

DATA PRODUCTS AVAILABLE FOR DOWNLOAD 
The following data products are available in the electronic version of this paper. 

1. Table [2] provides the set of 4714 radio quasars with visually-confirmed SDSS redshifts and radio morphological 
classifications. The main radio morphology classes are core, lobe, jet, triple, knotty. Example images are shown 
in Figure [U and the number of quasars in each class are listed in Table [H Triple quasars were sub-classified as 
core- dominated, lobe- dominated, or irregular, core quasars were sub-classified as resolved or unresolved. 

2. Table [3] provides a list of 317 double-lobed quasar candidates. These sources are quasars with SDSS spectra 
but with no radio detection within 2" of the optical position. They are identified by typically symmetric radio 
emission on either side of the core. Detailed selection criteria are discussed in Section [3731 
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Fig. 16. — The distribution of lobe-core-lobe opening angle and lobe-lobe flux density ratio for 519 triple-morphology quasars with 
well-determined lobe identifications (symbols). The dotted lines indicate the final selection criteria for this sample: an opening angle 
greater than 120° and a lobe-lobe magnitude difference of less than 2.5, corresponding to a ratio of < 10 in flux density. 



TABLE 9 

The 719 optically-faint triples in the SDSS 



SDSS position FIRST SDSS model magnitudes 



R.A. (J2000) Dec 


flux density [mjy] 


u 


g 


r 


i 


z 


2.020341 


-0.107590 


6.32 


22. Ti 


22.71 


21.97 


21.27 


20.62 


3.407161 


0.653116 


43.68 


25.13 


22.62 


21.86 


20.91 


20.67 


11.999740 


-8.890384 


2.23 


24.24 


23.59 


21.37 


20.11 


19.82 


25.386895 


1.011542 


3.49 


22.78 


21.69 


21.27 


21.13 
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26.709920 


0.628083 


2.91 


22.68 


22.50 


20.79 


19.78 


18.98 


28.002647 


-0.487744 


24.69 


23.08 


23.77 


23.10 


24.12 


20.95 


29.117733 


-9.749173 


28.06 


24.67 


22.03 


21.36 


21.02 


21.34 


29.719878 


1.025816 


59.42 


21.53 


20.81 


20.30 


19.96 


19.44 


32.471813 


-0.962489 


3.49 


23.68 


22.44 


22.23 


21.36 


20.97 


33.254807 


-0.304166 


52.80 


22.53 


22.59 


21.96 


21.77 


21.48 



Note. — This table is available in full in the electronic version of this paper. 

3. Table [9] lists 719 optically-faint triple-morphology quasars in the SDSS, with selection criteria described in 
the Appendix. These sources were identified purely on the basis of their radio emission, and have an SDSS 
photometric source within 2" of the radio core. These sources are not spectroscopically-confirmed quasars, as 
they are too faint for SDSS target ing, but many of them may qualiiy for spectroscopy in the SDSS III Baryon 
Oscillation Spectroscopic Survey ([Eisenstein et al.ll2011f) . 

4. Tables [S] and [5] provide median composite SDSS spectra for the radio-quiet quasars, for the entire radio quasars 
sample, for individual radio morphology classes (resolved core, unresolved core, triple, jet, lobe), and for triple 
quasars grouped by rest-frame R. 



Abazajian, K. et al. 2005, AJ, 129, 1755 



REFERENCES 

Abazajian, K. N. et al. 2009, ApJS, 182, 543 



22 



Kimball et al. 



Adclman-McCarthy, J. K. ct al. 2007, ApJS, 172, 634 

Antonucci, R. R. J.. & Ulvcstad, J. S. 1985, ApJ, 294, 158 

Baker, J. C. 1997, MNRAS, 286, 23 

Baker, J. C, & Hunstead, R. W. 1995, ApJ, 452, L95 

Baldwin, J. A. 1977, ApJ, 214, 679 

Barai, P., & Wiita, P. J. 2007, ApJ, 658, 217 

Barvainis, R., Lehar, J., Birkinshaw, M., Falcke, H., & Blundell, 

K. M. 2005, ApJ, 618, 108 
Becker, R. H., White, R. L., & Helfand, D. J. 1995, ApJ, 450, 559 
Bicknell, G. V. 1995, ApJS, 101, 29 

Blanton, M. R., Lin, H., Lupton, R. H., Maley, F. M., Young, 

N., Zehavi, I., & Loveday, J. 2003, AJ, 125, 2276 
Blundell, K. M., Rawlings, S., & Willott, C. J. 1999, AJ, 117, 677 
Brotherton, M. S., Tran, H. D., Becker, R. H., Gregg, M. D., 

Laurent-Muehleisen, S. A., & White, R. L. 2001, ApJ, 546, 775 
Browne, I. W. A., & Murphy, D. W. 1987, MNRAS, 226, 601 
Chiaberge, M., Capetti, A., & Cclotti, A. 1999, A&A, 349, 77 
CoUin-Souffrin, S. 1987, A&A, 179, 60 

Condon, J. J., Cotton, W. D., Crciscn, E. W., Yin, Q. F., Perley, 

R. A., Taylor, G. B., & Broderick, J. J. 1998, AJ, 115, 1693 
de Vrics, W. H., Becker, R. H., & White, R. L. 2006, AJ, 131, 

666 (dV06) 

de Vries, W. H., Becker, R. H., White, R. L., & Helfand, D. J. 

2004, AJ, 127, 2565 
Dietrich, M., Hamann, F., Shields, J. C, Constantin, A., 

Vcstcrgaard, M., ChafTcc, F., Foltz, C. B., & Junkkarincn, 

V. T. 2002, ApJ, 581, 912 
Dunlop, J. S., Taylor, G. L., Hughes, D. H., & Robson, E. I. 

1993, MNRAS, 264, 455 
Eiscnstcin. D. J. ct al. 2001, AJ, 122, 2267 
— . 2011. arXiv;1101.1529E 

Fanaroff, B. L., & Riley, J. M. 1974, MNRAS, 167, 31P 
Ferland, G. J., Peterson, B. M., Home, K., Welsh, W. F., & 

Nahar, S. N. 1992, ApJ, 387, 95 
Francis, P. J., Hewett, P. C, Foltz, C. B., Chaffee, F. H., 

Weymann, R. J., & Morris, S. L. 1991, ApJ, 373, 465 
Fukugita, M., Ichikawa, T., Gunn, J. E., Doi, M., ShimasaJcu, K., 

& Schneider, D. P. 1996, AJ, 111, 1748 
Ghisellini, G., Padovani, P., Celotti, A., & Maraschi, L. 1993, 

ApJ, 407, 65 
Gibson, R. R. et al. 2009, ApJ, 692, 758 
Gopal-Krishna, & Wiita, P. J. 2000, A&A, 363, 507 
— . 2001, A&A, 373, 100 

Gopal-Krishna, Wiita, P. J., & Dhurde, S. 2006, MNRAS, 369, 
1287 

Green, P. J. 1996, ApJ, 467, 61 
Gunn, J. E. et al. 1998, AJ, 116, 3040 
— . 2006, AJ, 131, 2332 

Hardcastle, M. J., & Worrall, D. M. 2000, MNRAS, 314, 359 
Harris, D. E., & Krawczynski, H. 2006, ARA&A, 44, 463 
Hogg, D. W., Baldry, I. K., Blanton, M. R., & Eisenstein, D. J. 

2002, ArXiv Astrophysics c-prints 
Hogg, D. W., Finkbcincr, D. P., Schlcgcl, D. J., & Gunn, J. E. 

2001, AJ, 122, 2129 
Hook, I. M., McMahon, R. G., Boyle, B. J., & Irwin, M. J. 1994, 

MNRAS, 268, 305 
Hornc, K., Korista, K. T., & Goad, M. R. 2003, MNRAS, 339, 

367 

Ivczic, Z. ct al. 2004, Astronomische Nachrichten, 325, 583 
— . 2002, AJ, 124, 2364 

Jackson, C. A., & Wall, J. V. 1999, MNRAS, 304, 160 
Jackson, N., & Browne, I. W. A. 1991, MNRAS, 250, 414 
Jackson, N., Browne, I. W. A., Murphy, D. W., & Saikia, D. J. 
1989, Nature, 338, 485 



Kapahi, V. K., & Saikia, D. J. 1982, Journal of Astrophysics and 

Astronomy. 3, 465 
Kelly, B. C. 2007, ApJ, 665, 1489 

Kharb, P., Lister, M. L., & Cooper, N. J. 2010, ApJ, 710, 764 
Kharb, P., & Shastri, P. 2004, A&A, 425, 825 
Kimball, A. E., & Ivezic, Z. 2008, AJ, 136, 684 {KI08) 
Korista, K., Baldwin, J., Ferland, G., & Vcrncr, D. 1997, ApJS, 
108, 401 

Ledlow, M. J., & Owen, F. N. 1996, AJ, 112, 9 

Lu, Y., Wang, T., Zhou, H., & Wu, J. 2007, AJ, 133, 1615 

Lupton, R. H., Ivezic, Z., Gunn, J. E., Knapp, G., Strauss, 

M. A., & Yasuda, N. 2002, in Proc. SPIE, 4836, 350 
Marscher, A. P., & Gear, W. K. 1985, ApJ, 298, 114 
McCarthy, P. J. 1993, ARA&A, 31, 639 
Morganti, R., Oosterloo, T. A., Reynolds, J. E., Tadhunter, 

C. N., & Migenes, V. 1997, MNRAS, 284, 541 
Morisawa, K., & Takahara, F. 1987, MNRAS, 228, 745 
Nenkova, M., Sirocky, M. M., Nikutta, R., Ivezic, Z., & Elitzur, 

M. 2008, ApJ, 685, 160 
Netzer, H. 1987, MNRAS, 225, 55 
Oke, J. B., & Gunn, J. E. 1983, ApJ, 266, 713 
Orr, M. J. L., & Browne, I. W. A. 1982, MNRAS, 200, 1067 
Osmcr, P. S., Porter, A. C, & Green, R. F. 1994, ApJ, 436, 678 
Osmer, P. S., & Shields, J. C. 1999, in Astronomical Society of 

the Pacific Conference Series, Vol. 162, Quasars and 

Cosmology, cd. G. Ferland & J. Baldwin, 235 
Padmanablian, N. ct al. 2008, ApJ, 674, 1217 
Peacock, J. A., Miller, L., & Longair, M. S. 1986, MNRAS, 218, 

265 

Pier, J. R., Munn, J. A., Hindslcy, R. B., Hcnnessy, G. S., Kent, 

S. M., Lupton, R. H., & Ivczic, Z. 2003, AJ, 125, 1559 
Rawlings, S., & Saunders, R. 1991, Nature, 349, 138 
Richards, G. T. ct al. 2002, AJ, 123, 2945 
— . 2003, AJ, 126, 1131 

Rys, S., & Machalski, J. 1990, A&A, 236, 15 
Schlcgcl, D. J., Finkbcincr, D. P., & Davis, M. 1998, ApJ, 500, 
525 

Schneider, D. P. ct al. 2007, AJ, 134, 102 
— . 2010, AJ, 139, 2360 
Scranton, R. et al. 2002, ApJ, 579, 48 
Shen, Y. et al. 2010, arXiv:1006.5178 
Shuder, J. M. 1981, ApJ, 244, 12 
Smith, J. A. et al. 2002, AJ, 123, 2121 
Strauss, M. A. et al. 2002, AJ, 124, 1810 

Tucker, D. L. et al. 2006, Astronomische Nachrichten, 327, 821 
Urry, C. M., & Padovani, P. 1995, PASP, 107, 803 
Vanden Berk, D. E. et al. 2001, AJ, 122, 549 
— . 2004, ApJ, 601, 692 

Weymann, R. J., Morris, S. L., Foltz, C. B., & Hewett, P. C. 

1991, ApJ, 373, 23 
Wills, B. J., & Brotherton, M. S. 1995, ApJ, 448, L81 
Wills, B. J., & Browne, I. W. A. 1986, ApJ, 302, 56 
Wu, J., Vanden Berk, D. E., Brandt, W. N., Schneider, D. P., 

Gibson, R. R., & Wu, J. 2009, ApJ, 702, 767 
Yee, H. K. C, & Okc, J. B. 1978, ApJ, 226, 753 
York, D. G. ct al. 2000, AJ, 120, 1579 
Zheng, W., Kriss, G. A., Tclfcr, R. C, Grimes, J. P., & 

Davidscn, A. F. 1997, ApJ, 475, 469 
Zirbel, E. L., & Baum, S. A. 1995, ApJ, 448, 521 



